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A weakly interacting massive particle (WIMP) weighing only a few GeV has been invoked as 
an explanation for the signal from the DAMA/LIBRA experiment. We show that the data from 
DAM A/LIBRA are now powerful enough to strongly constrain the properties of any putative WIMP. 
Accounting for the detailed recoil spectrum, a light WIMP with a Maxwellian velocity distribution 
and a spin-independent (SI) interaction cannot account for the data. Even neglecting the spec- 
trum, much of the parameter space is excluded by limits from the DAMA unmodulated signal at 
low energies. Significant modifications to the astrophysics or particle physics can open light mass 
windows. 



PACS numbers: 95.35. +d 

I. INTRODUCTION 

The DAMA/LIBRA Nal(Tl) scintillation experiment 
P] has used the annual modulation technique [2j [3] to 
search for dark matter (DM). They now find a modula- 
tion of over 8cr with a period and phase consistent with 
a DM signal. However, in some models of DM, it is not 
trivial to square this positive result with the null results 
from other direct detection experiments. 

Recent investigations [4] [5] [6] [7] , updating the discus- 
sion of |8j note that light DM, with mass of a few GeV 
might reconcile the DAMA/LIBRA data with constraints 
from other experiments, e.g., [5] [TO]. In addition, a DM 
candidate with mass ~ GeV is tantalizing - it might give 
insight into the ratio of the DM density to the density of 
baryons. 

In this note, we point out that the statistics of the 
DAMA/LIBRA data are now sufficiently powerful that 
an explanation of the DAMA/LIBRA data must now go 
beyond just fitting the overall modulation rate. Addi- 
tional self-consistency checks on the light WIMP scenario 
are now possible. We discuss two such checks. 

First, DAMA/LIBRA has now measured the modula- 
tion rate as a function of the observed recoil energy. This 
spectrum contains valuable information. Simple kine- 
matics indicate that the spectra of nuclei recoiling against 
a WIMP are sensitive to the mass of the WIMP. Thus, 
fitting the observed energy spectrum constrains the mass 
of the candidate WIMP particle. Another constraint can 
be derived by looking at the total (unmodulated) rate 
of observed events with low energy recoils. Some WIMP 
candidates will provide more events than the total num- 
ber of observed events at low energies, despite a presum- 
ably sizable background. 

These two constraints are powerful probes of the light 
WIMP region, and effectively exclude this interpretation 
if a Maxwellian velocity distribution with standard pa- 
rameters is assumed. Modifications to this assumption, 
in particular DM streams, can open up small regions of 
allowed parameter space. 



II. MODULATION RATES AND DAMA 

To calculate the detection rates, we follow the standard 
formalism reviewed in (TTJ [12] • The differential scattering 
rate is given by: 

dR ^ M N N TPx cr n F2 (f p Z+f n (A-Z)) 2 f f(y) 
dE R 2m x /4 e f 2 J Vm . n v 

(1) 

Here, Mm is the mass of the target nucleus, Nt is the 
number of target nuclei in a detector, and /j, ne is the 
reduced mass of the WIMP-nucleon system. We take 
the local DM density as p x = 0.3 GeV/cm 3 . F 2 is a 
nuclear form factor, while a n represents the cross section 
to scatter on a nucleon at zero momentum transfer. The 
relative coupling to protons and neutrons are given by 
f p and /„. We take these couplings equal. This choice 
does not have a large effect. We use the Helm form factor 
[TTJ [12]. Because the F 2 w 1 for the recoils in the low 
mass window, we do not expect our results to be sensitive 
to the detailed form of the form factor. 

The modulation signal arises from the seasonal differ- 
ences in the speed distribution, f(v). When the earth 
moves with the sun through the WIMP halo, the scat- 
tering rate above a given threshold is higher: the flux 
of WIMPs increases and the higher relative velocities fa- 
vor harder scatters. We take the WIMP halo to have 
a Maxwell-Boltzmann (MB) distribution with dispersion 
v = 220 km/s, cut-off by an escape velocity of 600 
km/s in the frame where the halo is isotropic. This 
value of the escape velocity is consistent with recent ob- 
servations from the RAVE survey [13] . which quotes a 
value of 498 km/s < v csc < 608 km/s gS]. The earth 
moves through this halo at approximately 240 km/s, so 
the highest velocity WIMPs in this boosted frame have 
velocities of approximately 840 km/s. In doing our nu- 
merical work, we take the earth velocity from [12j and 
the solar velocity relative to the rotation velocity (taken 
to be vo) from |14j . 

DAMA/LIBRA presents the modulation signal as a 
function of the observed energy, measured in keV elec- 
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FIG. 1: We show the region consistent with the DAMA/LIBRA modulation signal (at 68%, 90% and 99% CL), accounting 
for the detailed modulation energy spectrum between 2-6 keVee and the single bin modulation between 6-14 keVee (colored 
contours). Around this, we show the envelope (gray line) consistent with the DAMA/LIBRA modulation signal at 99% CL, 
fitting the total modulation rate from the two bins, 2-6 keVee and 6-14 keVee (i.e., no spectral information). Also shown are 
bounds calculated by considering data from the CDMS-II (Si) and XENON10 experiments. The limit curve labeled "DA MA- 
Total" arises from the self-consistency with the total (unmodulated) rates in the DAMA-LIBRA data, as described in Section III 
At right, we zoom in on the light mass region. When zooming in, we study variations about the minimum x 2 m the l° w mass 
region, as described in the text. 



tron equivalent (keVee). Conversion between true nu- 
clear recoil energy (keV) and observed energy (keVee) 
requires a quenching factor (QF^a =0.3, QFj= 0.09). 
We have included Gaussian energy smearing using the 
value: a(E)/E = 0.U8/VE + 9.1 x 10~ 3 (all energies in 
ke Vee) [T5] . It has been noted [5] that a low mass window 
consistent with other experiments, relies on the "channel- 
ing" effect [16 . When channeling occurs a recoiling nu- 
cleus interacts with the crystal only clcctromagnctically, 
efficiently transferring its energy to photoclcctrons. The 
result is that the full energy of a nuclear recoil is observed 
as scintillation light, rather than the quenched value. We 
include channeling following the parameterization of [7], 
also used in [5]. In this approximation, if channeling oc- 
curs in a given event, then all recoil energy is observed 
without quenching [26 . 

Fitting the DAMA modulation between 2-6 keVee and 
6-14 keVee as two bins, as done in j5j |8], we find a large 
region consistent at 99% CL, (gray unshaded contour in 
Fig. [TJ . There are two qualitatively distinct regions. One 
is centered roughly around 80 GeV, where the scattering 
dominantly proceeds off of the iodine nuclei. This re- 
gion is clearly excluded by other direct detection experi- 
ments. The second is at lower mass, extending down to a 
few GeV, where the scattering is dominantly off sodium 
nuclei. This DM is too light to cause recoils on iodine 
above detection threshold. Because a typical SI detection 
cross section scales with the square of atomic number, the 



cross section per nucleon for a WIMP that scatters off of 
sodium must be large. It is this low mass, high cross- 
section region that appeared consistent with all known 
direct detection experiments [27] . We now show that de- 
tailed spectral information modifies this region so that it 
is disfavored by current experimental results. 



III. DAMA VS. DAMA 

In Fig. [2] we show the observed modulation signal S m 
vs. observed energy for a WIMP with SI interaction at 
DAMA/LIBRA. Two of the plotted modulation spectra 
are for the masses that provide the best fits when the 
WIMP scatters dominantly on iodine (m x = 77 GeV), 
and on sodium (m x = 12 GeV). Our fit is to a nine 
bin x 2 from the eight .5 keVee bins between 2-6 keVee, 
where modulation is seen, and one bin from 6-14 keVee, 
where it is consistent with no modulation. The fit in the 
iodine scattering region is excellent (Xmm = 4-21 for 7 
dof, p = 0.76). The fit in the sodium region is worse: 
(Xmin — 7.54, p = 0.37). From Fig. [2J it is easy to see 
why this is so. At the best fit point in the low mass 
window, the DM is too light to have most of its recoils 
in the region of interest (2-6 keVee). The amount of 
modulation is still rising as one goes to lower energies 
- in conflict with the DAMA/LIBRA data. While one 
might think that this could be remedied by moving to a 
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FIG. 2: We show the modulation spectra for the best fit 
point where scattering off iodine dominates, m x = 77 GeV 
(dot-dashed orange), and three points where scattering off of 
sodium dominates. The best fit point off sodium is m x — 12 
GeV (solid red). We also show m x = 2 GeV (dashed green) 
and m x = 7 GeV (dotted blue). The points with error bars 
are the published DAM A/LIBRA data. 



higher mass, say 20 GeV, this approach does not succeed. 
As the mass moves above 12 GeV, a contribution coming 
from iodine scattering begins to move into the low end of 
the observed energy region, spoiling the fit. Also plotted 
in Fig. [2] are spectra for WIMP masses of 2 and 7 GeV 
for comparison. 

The 68%, 90%, and 99% CL (A X 2 < 2.3, 4.61, 9.21) 
contours consistent with our nine bin DAMA/LIBRA \ 2 
function are shown in Fig. [l] Both regions shrink dra- 
matically compared to the two bin \ 2 ■ I n the left panel, 
the A\ 2 is with respect to the global best fit point at 77 
GeV. In the right panel, we concentrate on the low mass 
region and have defined A\ 2 relative to the low mass best 
fit point of 12 GeV. Note this region is confined to masses 
above 10 GeV. This can be understood from examining 
the recoil spectra for the light WIMPs in Fig. [2] For a 
fixed overall modulation rate, sub- 10 GeV WIMPs pre- 
dict too little modulation above a couple of keVee: they 
simply do not have enough mass to cause recoils of this 
size. 

In Fig. [T| we have superimposed 90% limit contours 
from both the CDMS [H] and XENON [TO] experiments. 
We only show the CDMS contour relevant for low masses, 
corresponding to data taken with the silicon detectors in 
the Soudan mine, which have a 7 keV threshold. We have 
recalculated limits using the astrophysical parameters de- 
scribed here. For the XENON experiment, we account 
for the energy dependent efficiencies as described in |10j 
and to set the limit, apply the maximum gap method [T7] 
to the energy recoil range of 4.5-26.9 keV. We see that 
the light SI DAMA/LIBRA mass region is excluded once 
the modulation spectrum is taken into account. 

Another constraint can be imposed by looking at 
the total (unmodulated) counts at low energies at 
DAMA/LIBRA. For a WIMP, the number of recoil events 
increases with decreasing energy. Since DAMA/LIBRA 



cannot distinguish background events from signal events 
in this sample, it is clear that the predicted number of 
WIMP events should not excessively exceed their total 
number of counts in any bin. 

We require that the unmodulated rate in each bin from 
0.75 — 4 keVee not exceed the observed values within 
their 90% error [JJ. We show this constraint in Fig. [T] 
labeled "DAMA- Total." The allowed region lies below 
this curve. This constraint does not greatly impinge upon 
the allowed region from our nine bin \ 2 that accounts for 
the spectral details. Its constraint is most striking if one 
considers the allowed region of the two bin fit (see Fig. [IJ . 

Since the modulation in the low (2 — 2.5 keV) bin seems 
to be the most constraining, one can consider whether it 
is overly biasing our analysis. For instance, one could 
worry if DAMA/LIBRA were to restate the efficiency 
in the lowest bin, this might completely change our re- 
sults. We have explored such effects by various methods: 
tripling the error bar on the lowest bin, merging the en- 
tire range into a 2 — 3 keV bin, and discarding the lowest 
bin entirely. We find that only the last option (discard- 
ing the lowest bin) opens up a region of parameter space, 
with a point allowed with \ 2 = 9.14 for 6 dof (p = 0.17). 
This point also has a unmodulated rate that is close to 
saturating the observed rate. 



IV. VARIATIONS FROM ASTROPHYSICS AND 
PARTICLE PHYSICS 

Thus far, we assumed a MB halo and an elastic, SI in- 
teraction. Relaxing these assumptions could enlarge the 
region at light masses, so that modulation arises with an 
appropriate spectrum, consistent with other experiments. 

Let us begin by considering astrophysical modifications 
to the velocity distribution. Kinematics informs us of 
what modifications are needed. To scatter with nuclear 
recoil energy En, a WIMP must have a minimum veloc- 
ity (3 m i n = yj MnEr/2/j, 2 , where /1 is the reduced mass 
of the WIMP-nucleus (not nucleon) system. Consider 
the channeled possibility, for which the velocity require- 
ments are weakest: for scattering on sodium, with re- 
coil energy of 4.5 keV (the highest bin with significant 
modulation), one finds j3 m i n c ~ 1140,790,620 km/s for 
rn x = 2,3,4 GeV. If halo particle velocities approxi- 
mately follow a MB distribution, the most significant de- 
viations naturally occur for the highest velocities, where 
recent infall and streams may not have fully virialized. 
As such, the lightest particles are the most likely to have 
allowed regions opened by such deviations from a MB 
distribution. 

One modification to the halo is to include streams 
[TBI EH] ■ We investigated a wide range of streams, varying 
its velocity —1200 km/s < v str < 1200 km/s and disper- 
sion 10 km/s < <j str < 50 km/s. Our stream is such that 
for positive (negative) v s t r , the stream is directly against 
(with) the Earth's velocity as given in the sun's rest 
frame. Because we limit ourselves to small perturbations 



4 




FIG. 3: Similar to Fig.[T] but incorporating a stream of DM, 
with v 3tr = 900 km/s, and a 3 t T = 20 km/s. Insets: magnifi- 
cation of parameter space near 2 GeV and 4 GeV. The 2 GeV 
relies on both channeling and the stream, the 4 GeV region 
arises from unchanneled events from the stream. 



FIG. 4: Similar to Fig.[T] but including only inelastic scatter- 
ings, with S = 35 keV. 



enough to excite to the heavier state, i.e., 



on MB (we take the stream to have 3% of the density of 
the MB contribution), only streams well away from the 
bulk of the MB distribution make a difference. Shown in 
Fig. [3] for v str = 900 km/s, a str = 20 km/s, a channeled 
region near m x — 2 GeV appears with a marginal fit of 
X 2 = 10.3 and p — 0.17, as well as an unchanneled region 
near m x = 4 GeV with a fit of x 2 — 9.77 and p = 0.2. 
These low mass regions are allowed only for streams that 
are not bound to the galaxy, with v str > 800 km/s. As 
such, we view these modifications as significant from the 
typical MB distribution. Similar light regions for nega- 
tive v s tr are ruled out by the DAMA unmodulated limit. 

Direct detection experiments must go to very low 
thresholds to probe these theories. Comparisons within 
the context of MB halos will not be informative [2"B] . 
Rather, experiments should focus on probing a compara- 
ble range of velocity space. Probing the DAMA/LIBRA 
region of parameter space for nuclear recoils, requires 

Ethresh — 4 keV x where Mn is the target nu- 

cleus in question. For Si, a 3.3 keV threshold would be 
sufficient, while for Ge, and Xe, 1.25 keV and 0.7 keV 
thresholds arc required, respectively. Low thresholds for 
both Si and Ge are being explored |20j and given projec- 
tions from CoGeNT [2T], a recent low threshold germa- 
nium experiment, this low mass region could be probed 
even in the presence of non-MB distributions. 

A simple modification of the particle physics would be 
the inelastic DM scenario (iDM) [22], wherein scattering 
occurs via a transition to a slightly heavier state, m x * — 
m x = 5 ~ (3 2 m x . The WIMP must now be moving fast 



Such an inelastic WIMP at heavier masses (~ 100 GeV) 
has recently been shown to be consistent with DAMA's 
signal and the limits from other experiments [23 . 

Such scatterings have a recoil spectrum which can peak 
at 2-4 keVee, while being zero at low energies [55]. We 
have explored inelastic scenarios, and find regions with 
(5-30 keV, m x ~ 10 GeV and a n ~ KT 39 cm 2 that 
are consistent with other experiments. XENON10 is the 
most constraining. We show the allowed region for 5 — 
35 keV in Fig. gj 

Building an inelastic model with such a large cr n is 
difficult. In the prototypical inelastic model, scattering 
arises from Z exchange. A light, vector- like Dirac neu- 
trino would scatter off sodium with a n « 1.5 x 10 -39 cm 2 . 
But since the number of neutrinos is limited by LEP to 
be n v = 2.9840±0.0082 ;24], the Z-DM coupling must be 
suppressed, for instance by mixing with a singlet state. 
This bounds the cross section cr„ ^ 1.3 x 10~ 41 cm 2 . 
Thus, there is significant tension between a viable model 
of light inelastic DM, and the overall signal size. 

Nonetheless, modifications to the astrophysics and par- 
ticle physics can re-open some windows in the light 
WIMP regime, with astrophysics likely opening m x — 
2 — 5 GeV regions, and inelasticity opening m x — 10 GeV 
windows. Both of these warrant further investigation. 
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V. CONCLUSIONS 

The recent update from the DAMA/LIBRA collabo- 
ration has added significance to the modulated signa- 
ture. Detailed spectral data should now be considered 
for any model proposed to explain the signal. As it ap- 
pears that modulation extends up to at least 5 keVee, 
the spectrum places very strong constraints on the prop- 
erties of a standard WIMP in a MB velocity distribution. 
Additionally, because the signal from a standard WIMP 
rises very rapidly at low recoil energies, one must insure 
that the unmodulated signal predicted does not exceed 
the observed levels at DAMA/LIBRA. 

Fitting to the spectral data, one finds two regions for 
a standard WIMP in a Maxwellian halo which explain 
the spectrum. The high mass region conflicts with other 
direct detection experiments and the unmodulated low- 
energy rate at DAMA/LIBRA. The light mass region 
gives a poorer fit to the spectrum, but is still acceptable 
and is largely consistent with limits from the DAMA un- 
modulated signal. However, 90% CL limits from CDMS- 



Si and XENON-10 exclude this region. 

Modifications to astrophysics or particle physics can 
open light mass windows. Streams can open regions in 
the m x ^2 — 5 GeV range, although for very particular 
choices of parameters. Inelasticity can broaden areas at 
higher masses, but faces model-building challenges. Fu- 
ture direct detection experiments should push to explore 
both of these regions. 
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